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Abstract 


Anode two-phase flow has an important influence on DMFC performance and methanol crossover. In order to elucidate two-phase flow influence 
on anode performance, in this work, anode polarization is investigated combining experimental and modelling approach. A systematic experimental 
analysis of operating conditions influence on anode polarization is presented. Hysteresis due to operating condition is observed; experimental results 
suggest that it arises from methanol accumulation and has to be considered in evaluating DMFC performances and measurements reproducibility. 
A model of DMFC anode polarization is presented and utilised as tool to investigate anode two-phase flow. The proposed analysis permits one to 
produce a confident interpretation of the main involved phenomena. In particular, it confirms that methanol electro-oxidation kinetics is weakly 
dependent on methanol concentration and that methanol transport in gas phase produces an important contribution in anode feeding. Moreover, it 
emphasises the possibility to optimise anode flow rate in order to improve DMFC performance and reduce methanol crossover. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Hydrogen polymer electrolyte membrane fuel cell (PEMFC) 
is considered a highly promising technology, especially for 
micro power generation and vehicular applications, due to its 
important attributes: low temperature and low pressure oper- 
ation, no liquid electrolyte, high power density. The direct 
methanol fuel cell (DMFC) technology is a further develop- 
ment of PEMFC, particularly promising for portable electronics 
and vehicular applications, due to the advantage of a liquid 
fuel. Its main drawbacks are efficiency and power density lower 
than PEMFC. These topics are intensively studied [1,2] and 
are mainly attributed to methanol permeation through the poly- 
mer membrane and slow electrochemical methanol oxidation. 
The influence of anode two-phase flow is generally less inves- 
tigated, despite its effect on DMFC performance and methanol 
crossover flux is very important. Experimental studies indicate 
degradation of cell performance under intensive generation of 
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bubbles in the flow field [3-5]. Different works investigated 
operating conditions influence on anode polarization [7—9], 
but generally measurement uncertainty and reproducibility are 
not evaluated. Due to the complexity and the time needed 
for a whole experimental investigation numerical models have 
been adopted by some authors to investigate DMFC behaviour. 
But also recent state-of-the-art numerical models of DMFC 
present difficulties in rigorous description of the two-phase flow 
[10,11]. Kulikovsky developed a DMFC model [12], describ- 
ing the two-phase flow in the anode channel, with an adaptive 
parameter, showing that under moderate methanol and oxygen 
stoichiometry the bubbles reduce limiting current density of the 
cell considerably. Ge presented a three-dimensional two-phase 
flow model [13], assuming only carbon dioxide presence in 
anode gas phase, and showed the influence of anode two-phase 
regime on methanol crossover. Wang presented a complete two- 
dimensional two-phase model, solved with a CFD technique, 
validated on a limited set of experimental measures [14]. 

This work aims to investigated anode polarization, in particu- 
lar to produce a better understanding of two-phase flow influence 
on anode performance. A combined experimental and modelling 
approach is proposed. A systematic experimental analysis is 
reported, measurements are traceable to the international ref- 
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Nomenclature 

Cref reference concentration (mol cm~?) 

C species concentration in channel (mol cm7?) 

C time-averaged concentration in channel 
(mol cm~?) 

Ci species concentration in catalyst layer (mol cm~?) 

Dy effective diffusivity in diffusion layer (cm? s7!) 

Dm effective diffusivity in membrane (cm? s7!) 

Eo ideal potential difference (V) 

F Faraday constant (C mol`!) 

h channel height and width (cm) 

i local current density (A cm~?) 

ix exchange current density (A cm7?) 

k Tafel constant (A cm7?) 

Ky Henry constant (mol J7 9) 

lb diffusion layer thickness (cm) 

A catalyst layer thickness (cm) 

lm membrane thickness (cm) 

L channel length (cm) 

M molecular weight (g mol~!) 

m inlet mass flow rate (g min—!) 

Neross crossover flux (mol cm~? s7!) 

nd drag coefficient 

Psat saturation pressure (MPa) 

P channel pressure (Pa) 

R universal gas constant (J mol`! K7!) 

Re contact resistance (Q cm?) 

T fuel cell temperature (K) 

v local velocity in channel (cm s7!) 

Veell cell potential difference (V) 

x x coordinate (cm) 

y y coordinate (cm) 

z z coordinate (cm) 

Greek symbol 

a Tafel transport coefficient 

€ volumetric void fraction 

n polarization (V) 

y reaction order 

p compound density (g cm7?) 

Om membrane conductibility (Q7! cm7!) 

Superscript 

a relative to anode 

G relative to gas phase 

L relative to liquid phase 

Subscript 

CO2 relative to carbon dioxide 

H20 relative to water 

in relative to channel inlet 

Met relative to methanol 

out relative to channel outlet 


erence standard and are characterised in terms of uncertainty. 
This complete experimental characterisation permits to validate 
part of the model published in [15], analysing accurately model 
residuals and comparing them with experimental uncertainty. 

Both the experimental analysis and the accurate model val- 
idation permit to produce a solid interpretation of the main 
phenomena involved in DMFC operation, contributing to opti- 
mise operating conditions and components. 


2. Experimental equipment 


The experimental setup is composed of a DMFC hardware, 
having a maximum cross-sectional area of 25 cmĉ?and a unique 
membrane electrode assembly (MEA) with an effective area of 
12.25 + 0.5 cm?. The MEA is contained between two graphite 
blocks in which flow distributors are cut out (single serpentine 
channel, square section: depth 0.8mm, width 0.8 mm, length 
1500 mm). The cell is held together with two stainless steel plates 
using eight retaining bolts, which were closed applying a torque 
of 12+0.5Nm. 

To accommodate a calibrated thermocouple (uncertainty 
0.05 K), connected both with a temperature controller and an 
acquisition system, a slot is present in the cathode steel plate. 
Two electrical heaters are placed within the steel plates, one 
for each electrode, to control the cell temperature. The graphite 
blocks present electrical contacts with gold plates, connected to 
a high accuracy power supply (current uncertainty 0.1% of the 
effective value + 0.012 A; voltage uncertainty 0.05% + 0.008 V). 
The anode solution is fed by a peristaltic pump with a resolution 
of 1 rpm and a speed uncertainty of 0.5%; anode mass flow rate 
is measured with an uncertainty of 1%. The methanol solutions 
were produced once, at the beginning of the experimentation, 
mixing bidistilled water and methanol (grade 99.5%mass) and 
measuring the solution mass fractions (uncertainty 0.05 %mass) 
with a calibrated balance (uncertainty 0.1 g). Hydrogen, pro- 
duced by a water electrolyzer (grade 99.999 % mole), is supplied 
to cathode side, as a reference electrode; its flow rate, constant 
at 0.2 N dm? min! is controlled and measured by a calibrated 
flowmeter (uncertainty 0.002 N dm? min~!). 

The MEA membrane is GEFC-117, anode catalyst layer 
(thickness 10 um) presents a metal loading of 4mg cm? 
(Pt:Ru=2 on weight basis), cathode catalyst layer (thickness 
10 um) presents a metal loading of 4mgcm? (Pt), anode and 
cathode diffusion layers are carbon paper (thickness 100 wm), 
cell gaskets consist of a Mylar layer (thickness 100 um), in 
contact with the membrane, and a fibreglass layer (thickness 
150 wm), covered by PTFE. The MEA was activated through 
a conditioning procedure, that consists of operating periods 
of 30min at constant temperature (343 K) and load (0.6 V), 
with feeding of hydrogen (0.2 + 0.002 N dm? min7!) and air 
(2+0.02N dm? min~!) saturated by water at 353 K and ambi- 
ent pressure, for a total duration of 18 h subdivided equally in 3 
days. 

The described equipment permits thus to determine anode 
polarization, because cathode is considered a reference elec- 
trode, in particular a pseudo-DHE (dynamic hydrogen electrode) 
[6], whose potential is nearly constant and negligible in com- 


374 A. Casalegno, R. Marchesi / Journal of Power Sources 175 (2008) 372-382 


parison to anode, when determining anode polarization curves. 
Therefore, measured voltage variations are mainly attributable 
to anode and membrane overpotentials, as described in [7], and 
the following polarization curves are reported versus reversible 
hydrogen electrode. 

An acquisition procedure was developed to determine anode 
polarization characteristics for defined operating conditions and 
to permit measurement reproducibility evaluation. It includes 
two steps: initial transitory and characteristic acquisition. The 
initial transitory consists in an acquisition period, performed in 
the investigated operating conditions, fixing the voltage at 0.5 V, 
measuring voltage and current at 1 Hz frequency for 15 min; it 
is preceded by a preliminary not acquired transitory period of 
15 min. The characteristic acquisition is composed of 22 single 
point measurements referred to different voltages (from 0.25 
to 0.75 V, stepped by 0.05 V), collected following two one-way 
curves respectively increasing and decreasing voltage. Each sin- 
gle point acquisition is performed at constant voltage, measuring 
voltage and current at 1 Hz frequency for 400 s. Reproducibility 
procedure is composed of an initial transitory and two single 
point acquisitions, performed at 0.4 and 0.6 V. 

Data obtained by a single point acquisition (two series of 400 
values, voltage and current) are elaborated in two steps: transi- 
tory elimination and outliers elimination. The first 60 values are 
discarded, because resulting from the transitory caused by volt- 
age variation. A robust method is used to individuate outliers. It 
eliminates values not included in the interval median + 3 times 
standard deviation, estimated through median absolute devia- 
tion. Finally, the last 180 points are selected as significant, then 
mean values of voltage and current are calculated [16,17]. 

Measurement uncertainty of current is evaluated according to 
[16,17]. The contribution of each controlled parameter is eval- 
uated considering the uncertainty of each measuring device and 
estimating experimentally its influence on current. The final cur- 
rent uncertainty results in the geometric sum of the uncertainties 
produced by the different controlled parameters. The single point 
uncertainty, u;, (179 dof, 95% population) for current density i 
from 0 to 1.1 A cm7?, is estimated equal to: 


uj = —5.8 x 10-777 + 6.8 x 10-71 +6 x 107° (1) 
3. Analysis of experimental results 


The following work refers to the experimental results of 
eight anode characteristic curves, whose operating conditions 
are reported in Table 1. The aim of the analysis is the investiga- 
tion of operating conditions influence on anode polarization. 


3.1. Hysteresis 


Comparing two one-way characteristics, obtained increas- 
ing and decreasing voltage as previously described, a hysteresis 
phenomenon is observed; two examples are shown in Fig. |. For 
high voltages the difference between the two average current 
densities has a magnitude order of 1072A cm~?, similar to the 
typical uncertainty of single point current measurement. In some 
cases, the differences are statistically significant, as confirmed 


Table 1 
Controlled parameters of experimental sessions 


Session Methanol concentration Cell temperature Anode flow rate 
(wt.%) (uncertainty (K) (uncertainty (g/min) 
0.07 wt.%) 0.05 K) (uncertainty 1%) 

1 3.25 332.9 0.47 

2 3.25 333.1 1.08 

3 3.25 353.6 0.47 

4 3.25 353.3 1.08 

5 6.5 332.7. 0.47 

6 6.5 332.8 1.08 

7 6.5 352.7 0.47 

8 6.5 352.7 1.08 


by the analysis of variance (ANOVA), and the two mean values 
have to be considered different. The differences are significant 
exclusively at high current density, in particular where fuel util- 
isation is above 80%. It is interesting to observe that during the 
backward curve, decreasing voltage, the current is lower than 
during the forward curve. The hysteresis may be caused by the 
following three factors: methanol, CO2 or water accumulation. 
In the complete set of experimental data, the hysteresis results 
statistically significant exclusively when the fed solution has a 
concentration of 3.25 wt.%. 

With a solution concentration equal to 6.5 wt.% instead the 
hysteresis is not significant, as evident in Fig. 2. The curves at 
6.5 wt.% reach higher current densities, thus a more intense CO2 
production and concentration, implying that hysteresis is not due 
to CO accumulation, that in principle could reduced active fuel 
cell area. Moreover, at higher current density void fraction at 
anode increases, reducing water availability, allowing to exclude 
that hysteresis is caused by local membrane dehydration. 

The previous considerations seem to individuate the cause 
of hysteresis in some sort of long-term methanol accumula- 
tion, coherently with [18]: during forward curve accumulated 
methanol is available to be consumed, on the contrary during 
backward curve there is methanol lack. This conclusion does 
not pretend to be exhaustive, further experimental analyses are 
needed to produce a more rigorous quantitative demonstration. 
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Fig. 1. Hysteresis in polarization curves. Fuel cell temperature: 333 K; methanol 
concentration: 3.25 wt.%. 
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Fig. 2. Hysteresis in polarization curves. Fuel cell temperature: 333 K; methanol 
concentration: 6.5 wt.%. 


Consequently hysteresis depends on operating condition history 
and have to be considered in evaluating the uncertainty of the 
mean values between forward and backward curves. Prudently 
final current density uncertainty, Eq. (1), is expanded by a factor 
equal to 1.7, compared to single point uncertainty; in this way 
ANOVA confirms that hysteresis effect is no more statistically 
significant. 


3.2. Reproducibility 


Reproducibility analysis is necessary to characterise properly 
the measures and to verify their reliability. Five reproducibil- 
ity measures were carried out each in a different day for every 
polarization curve reported in this work; in Fig. 3 an example is 
reported, indicating current density measures in different days, 
their average value and uncertainty. 

Current density variation between different days has gen- 
erally a magnitude order of 1073 A cm~?, inferior than single 
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Fig. 3. Reproducibility measures in five different days. Fuel cell temperature: 


353 K; methanol concentration: 6.5 wt.%; anode flow rate: 1.08 g min™!. 
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Fig. 4. Polarization curve. Methanol concentration: 3.25 wt.%; varying anode 
flow rate and fuel cell temperature. 


point uncertainty. In few cases, around 10% of the analysed data, 
ANOVA reveals significant differences between days, consider- 
ing single point uncertainty, not expanded by hysteresis factor. 
Taking into account the prudently adopted hysteresis factor, no 
differences between days are significant and reproducibility is 
demonstrated. 


3.3. Operating conditions influence 


In Figs. 4 and 5, four polarization curves at fixed methanol 
concentration are compared. The expected influence of temper- 
ature and anode flow rate is observed [7—9,19,20]: increasing 
temperature, methanol oxidation kinetics is enhanced, due to a 
reduction of activation losses; increasing anode flow rate, limit- 
ing current is higher and concentration losses reduce. Moreover, 
it is interesting to evidence that at 353 K the limiting currents 
result higher at the same anode flow rate and methanol con- 
centration, thus fuel utilisation increases. This behaviour can be 
explained considering that at high temperature methanol evap- 
oration and transport through gas phase is enhanced, increasing 
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Fig. 5. Polarization curve. Methanol concentration: 6.5 wt.%; varying anode 
flow rate and fuel cell temperature. 
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Fig. 6. Polarization curve. Fuel cell temperature: 353 K; varying methanol con- 
centration and anode flow rate. 


methanol availability, despite its crossover increases due to an 
increase in membrane permeability. 

Varying anode flow rate and methanol concentration at a fixed 
temperature, Fig. 6, it is evident that methanol electro-oxidation 
kinetics is weakly dependent on methanol concentration. 
Despite important variations in methanol concentration and 
anode flow rate, the polarization curves are very similar at 
low current, confirming that activation losses do not depend 
strongly on methanol concentration [8,9,19]; the main observed 
effects are due to concentration losses at high current density 
and methanol utilisation. 

In Fig. 7, anode performances at constant effective methanol 
flow rate, defined as the product of anode flow rate times 
methanol concentration, are compared. At a constant tempera- 
ture of 353 K the curves are very similar, implying that methanol 
availabilities in the two feeding configuration are comparable. 
At 333 K the behaviour is different: a higher methanol con- 
centration permits to reach higher limiting current. This can 
be interpreted considering that at lower temperatures methanol 
availability in gas phase is reduced, confirming the important 
role of gas phase methanol transport in feeding the electrode. 


0.8 T T r T T T 
0.7 + | 
= 0.6 + | 
ei 
fa 
& 0.55f 1 
© 
> 
0.4 + A 
—O— 333 K, 6.5%w, 0.47 gmin 
—t-— 333 K, 3.25%w, 1.08 gmin“! 
~~ —+- 353 K, 3.25%w, 1.08 gmin“ 
—*— 353 K, 6.5%w, 0.47 gmin™' 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 
Current Density [A-cm?] 


Fig. 7. Polarization curve. At similar effective methanol flow rate. 
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Fig. 8. Efficiency. Fuel cell temperature: 353 K; varying methanol concentration 
and anode flow rate. 


The fuel cell efficiency of methanol electrolysis can be 
defined, without considering eventual methanol recirculation, 
as the ratio between the effective produced power, equal to pro- 
duced hydrogen power minus electrical power consumption, and 
the power related to inlet methanol flow rate: 


(Eo — Veen) 


=r (2) 
nC? o LHV met 


The maximum efficiency reached in the set of experimental 
data is around 45%, Fig. 8, without considering eventual fuel 
recirculation. Moreover, at similar methanol flow rate the effi- 
ciency results higher at a higher concentration and lower anode 
flow rate, confirming that the losses typically observed at high 
concentration in fuel cell operation are only due to cathode 
polarization, caused by methanol crossover. 


4. Model validation 
4.1. Model description 


4.1.1. Assumptions 

The model does not pretend to be a rigorous and exhaustive 
description of complex DMFC behaviour, but aims to reproduce 
the main involved phenomena and consequently to be a tool for 
DMFC understanding and optimization. The model is based on 
the following assumptions, already reported in [10-14,23-25]: 


(1) The fuel cell hardware to be modelled has single serpen- 
tine square section distributors. The modelling domain is 
simplified as: a straight channel, anode distributor, and the 
MEA fed by it, Fig. 9. Coherently with the hardware anode 
gas channel height and width are equal, h*, and MEA 
width is two times h*. The MEA is composed by three 
sub domains: anode diffusion layers, catalyst layers and 
the membrane. The overall MEA area is assumed active. 
No reactions occur at cathode side. 

(2) The modelling domain is assumed isothermal. 

(3) The model neglects 2D and 3D effects; only convective 
transport is considered in the channel, along x direction; 
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diffusionlayer 
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diffusionlayer 
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= Cathode channel = 


Fig. 9. Fuel cell scheme. 


only diffusive transport takes place in y direction within 
diffusion layer and membrane; concentrations in catalyst 
layer are assumed constant. 

(4) Three compounds, methanol, water and carbon dioxide, 
are considered in two phases at anode side. 

(5) The anode two phase flow is simplified as a plug flow, 
where plugs and bubbles have the same velocity, constant 
within yz planes, as are species concentrations. Due to 
alternating plugs and bubbles single species concentrations 
in a certain volume are rapidly varying as time function, 
from the concentration in the liquid plug to the concentra- 
tion in the bubble. Time-averaged local concentration of 
methanol, water and carbon dioxide, indicated as C, are 
introduced to describe this behaviour at anode side. 

(6) Gases, liquids and mixtures have an ideal behaviour. At 
anode side the concentration of the species, present in gas 
and liquid phases, are in equilibrium. Diffusion and cata- 
lyst layer porous structure do not influence single species 
gas-liquid equilibrium and void fraction. 

(7) Anode pressure is assumed constant. 

(8) The membrane is considered perfectly hydrated, due to 
high water stoichiometry at anode side. 

(9) Methanol and water crossover are due to diffusion through 
the membrane and electro-osmosis. Methanol and water 
diffusion is promoted by their concentrations in liquid 
phase, as some experimental analysis suggest [7,21]. Mem- 
brane permeability to gaseous compounds is neglected, in 
particular carbon dioxide permeation through the mem- 
brane is negligible [22]. 

(10) Electrochemical reaction kinetics obey Tafel equation. 


4.1.2. Anode channel 

Taking into account the presented assumptions the species 
equations for methanol, carbon dioxide and water assume the 
following forms: 


hè d(v*Cuet) i 
2 ax =s 6F Nevo GB) 


h? alv’ Co,) i 


= 4 
2 ox 6F 4) 
AP Ciao) __t__yigo (5) 
2 ax = 6F cross 


where v* is the local velocity in x direction. Locally the con- 
vective flow of each species has a gradient equal to the sum of 
fluxes in y direction, caused by the anodic electrochemical reac- 
tion, in function of current density i, and by the crossover flux 
through the membrane. The factor h®/2 is due to assumption (1) 
and determines that MEA active area is equal to 2h*L. Bound- 
ary conditions are related to species inlet flows, their values are 
reported in Table 3. 

Time-averaged local concentration of methanol, water and 
carbon dioxide are defined as the sum of gas phase concen- 
tration, CS, times volumetric gas fraction € and liquid phase 
concentration, CH, times volumetric liquid fraction 1 — e€: 


CMet = éa F dd = Cris (6) 
Céo, = eCo, + 1 — OCG, (7) 
Cho = CHo + (1 — OCH (8) 


Single species concentrations in gas and liquid phase are con- 
sidered in equilibrium, assumption (6), described by Henry’s 
law for methanol and carbon dioxide. Moreover, gas phase is 
considered completely saturated by water: 


Chet = KU MetC Met RT (9) 
ch. =K CSa, RT (10) 
CO, H,CO2 CCo, 

H20 

Ds 
Cho = “RT (1) 


Ideal gas mixture is assumed for gas phase, (6), thus, defining 
P* as constant anode pressure, (7), Dalton equation relates gas 
phase concentrations: 

a 


P 
Chet + Cho + Ceo, = pp (12) 


Ideal liquid mixture implies that total mass fraction in liquid 
phase is equal to 1: 


Moet Mpo Mco 
L e a,L 2 a,L 2 — 
CMet pA + Cho T Cn e = 1 (13) 


pĂ is total density of liquid phase. It can be reasonably approx- 
imated equal to pure water density, Pio» considering that pure 
compounds densities are not very different from water value and 
compounds concentrations are limited to low values. 


4.1.3. Anode diffusion layer 

The sum of methanol fluxes in x direction, (3), is equal to 
the diffusive flow through the diffusion layer. The latter can be 
expressed as the sum of diffusive fluxes in both gas and liq- 
uid phases, assuming that void fraction is constant in y, within 
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diffusion and catalyst layers, (6): 


i Met De L L 
6F ae Nexoss = A (CMet Ci ) — €) 
b 
pèc 
G G 
t-g (Cha — CDE (14) 


where De and Dp are effective methanol diffusivity in liq- 
uid and gas phases within the diffusion layer, ce and Cc are 
methanol concentrations in the catalyst layer. 


4.1.4. Anode catalyst layer 

Time-averaged local concentration of methanol and gas- 
liquid equilibrium in the catalyst layer are defined similarly as 
in (6) and (9): 


Či = «cf +d=9ct (15) 
Ch = KymeCORT (16) 


Anode polarization is calculated assuming Tafel kinetic, (10), 
of order y*, depending on time-averaged local concentration of 
methanol in the catalyst layer: 


foe T S (17) 
= atk ka 
cata 
e=a($) (18) 
ref 


4.1.5. Membrane 

Methanol and water crossover are due to liquid phase diffu- 
sion and electro-osmosis, assumption (9). Water and methanol 
concentrations at cathode side are negligible, in comparison to 
anode concentrations. 

Electro-osmotic flux is proportional to current density, water 
drag coefficient of proton, na, and molar fraction of the involved 
species. Considering the abundant presence of water in the anode 
channel and the limited methanol concentration, the gradient of 
water concentration through the diffusion layer can be neglected. 
Thus, for simplicity, water concentration in the catalyst layer is 
assumed equal to anode channel one. Moreover, for the purpose 
of electro-osmotic flux calculation, water molar fraction in the 
catalyst layer is assumed equal to 1. 


pMet i cl 
Novos = Cr + aE ar (19) 
m CHo 
DIO i 
Neos = 7 Cho + "az (20) 


Proton flow through the polymeric membrane produces a polar- 
ization, that obeys to the Ohm’s Law; thus, considering the 
membrane perfectly hydrated and its conductivity, om, constant 
in membrane domain: 


lm 
Nohm = I (=) (21) 


Table 2 

Geometrical parameters 

h? (cm) 0.08 
h! (cm) 0.08 
L (cm) 76.6 

l (cm) 0.01 
I, (cm) 0.01 
i} (cm) 0.001 
K (cm) 0.001 
Im (cm) 0.018 
Table 3 

Experimental conditions 

P* (Pa) 101325 
P® (Pa) 101325 
T(K) 333, 353 
m? (g min!) 0.47, 1.08 
m° (dm? min™!) 0.2 

C? et (Wt. 7%) 3.25, 6.5 
Veell (V) 0.25-0.75 


4.1.6. Cell potential difference 

Potential difference across the cell is constant along z direc- 
tion, as effect of electrodes equipotentiality, and is equal to the 
sum of anode, ohmic polarization and contact resistance voltage 
drop!: 


Veell = n? + Nohm + Rci (22) 


4.1.7. Numerical solution 

The system, composed by the equations from (3) to (22), 
presents three differential equations, 20 algebraic equations, 20 
variables. This DAE system can be solved numerically, apply- 
ing the appropriate initial conditions, regarding inlet flows, fuel 
cell temperature and potential difference. Total and mean elec- 
trical current can be calculated by integration over MEA area. 
These quantities permit to build polarization curves, commonly 
used to characterise direct methanol fuel cells. For this work 
Matlab® environment was used to solve the DAE system. 


4.2. Validation 


The quantities present in the model have physico-chemical 
meaning, there are no pure adaptive parameters. These quan- 
tities are divided in: geometrical parameters, defined by the 
experimental hardware, shown in Table 2; operating conditions, 
defined by experimental analysis, presented in Table 3; assumed 
parameters, provided by literature, Table 4; fitting parameters, 
Table 5. 

The calibration procedure consists in minimisation of the 
residuals between model estimation and experimental results 
in the complete investigated range of operating conditions. The 
utilised experimental data are composed of 88 current density 


' Contact resistance and cable voltage drop are measured and reported in 
Table 4. 
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Table 4 
Assumed parameters 
a? Kı Calibration 
Ceg (molem 3) 1 x 1073 [14,25] 
DL (cm? s7!) 0.7 x 10714163 — (999.778/T) x 104 [14] 
Dy G (cm? s7!) Ks Calibration 
Dm (em? s~!) 1.28 x 1076 el (2416/303)—(2416/T)) [26] 
Eo (V) 1.2 
F (Cmol"!) 96495 
ys K4 Calibration 
i? (A cm~?) Ky e(£3/RX0/353)-0/T)) Calibration 
Ku Me (mol rh 2.2 @5200((1/T)—1/298)) [27] 
Ky,co, (mol J) 3.5 e2400((1/T)—(1/298)) [27] 
Na 2.9 @1029((1/333)—(1/T)) [23] 
Psat (MPa) —7.9((373/T) — 1) — 1.38 x 10-7(10!! 340-7373) — 1) + 8.13 x 1073(107349673/D)-D — 1) + [27] 

5.03 log(373/T) + log 1013 
R (Jmol~! K7!) 8.314 
R: (Q cm?) 7.5 x 107? Measure 
p“ (gem) = Pho 
peg (gcm?) (—0.0028(T — 273)? — 0.1757(T — 273) + 1003.8) x 107° [27] 
om (Q7! cm7!) 0.073 ef !268((1/298)-1/7))] ae 


measures, coming from eight polarization curves, previously 
presented, at different temperature, inlet methanol concentra- 
tion and inlet anode flow rate. The model is not intended to 
describe unsteady state behaviour, including hysteresis, thus its 
estimations are compared to mean data, characterised by final 
current density uncertainty, expanded by the hysteresis factor. 
Considering that the potential difference range interesting for 
applications is from 0.25 to 0.5 V, the calibration procedure 
was developed to produce a more accurate fitting in this range, 
where two-phase flow and methanol gas transport effects are still 
present. Moreover, at higher voltages complex two-phase fluid 
dynamics effects, not described by the model, may occur. 

The fitting parameters, reported in Table 5, characterise anode 
kinetics and effective methanol diffusivity in gas phase through 
diffusion layer. The fitting parameters, Kı, K2 and K3, char- 
acterise anode kinetics, composed of quantities affected by 
high uncertainty, demonstrated by a high variability in literature 
[12,14,23,24]. The values resulting from calibration are reason- 
able and very similar to literature, in particular K3 with [19,20]. 
The order of the overall anode reaction, K4, is another uncertain 
parameter, in literature varies from 0 [14] to 0.5 [12] to 1 [11]. 
The value obtained by calibration, 0.02, is consistent with the 
experimental observation of a very weak current dependence on 
methanol concentration in the activation region, Fig. 6. 

Considering the contribution of CO2 counter diffusion, vari- 
ation in void fraction and in methanol state equilibrium, the 
effective methanol diffusivity in gas phase through diffusion 
layer, K5, is a considerably uncertain parameter. The fitting 


Table 5 

Fitting parameters 

Kı 0.575 

Ka (A cm~?) 2.38 x 1074 
K3 (Jmol!) 73167 

K4 0.02 

Ks (cm? s7!) 3 x 1073 


values result 40 times lower than that assumed in [14]. This dif- 
ference is acceptable, considering the contribution and influence 
of the mentioned phenomena, moreover, the value is consis- 
tent with that obtained with the presented model and different 
experimental analysis, in [15]. In spite of a lower value of effec- 
tive methanol diffusivity than expected, the model confirms that 
methanol diffusion in gas phase has a considerable role in trans- 
port through the diffusion layer. 

The residuals produced in the investigated interval, 0.25-0.5 
V, between model estimation and experimental results have rea- 
sonably normal distribution and an absolute value inferior than 
0.03A cm~?, an example is reported in Figs. 10 and 11. 

An accurate residual analysis has been effectuated, evaluating 
them in comparison with experimental uncertainty. In Fig. 12, 
the residuals are reported in function of current density, with 
relative experimental uncertainty. Model accuracy is evaluated 
quantitatively by F-test, as reported in [28]. It verifies if the 
model is sufficiently accurate, in comparison to experimental 
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Fig. 10. Polarization curve experimental (points) and modelling (lines). Fuel 
cell temperature: 333 K; varying methanol concentration and anode flow rate. 
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Fig. 11. Polarization curve experimental (points) and modelling (lines). Fuel 
cell temperature: 353 K; varying methanol concentration and anode flow rate. 


uncertainty, to permit to consider the latter as confident interval 
of model estimation. 

Experimental measures from 0.25 to 0.5 V satisfy F- 
test (95%), thus in this range the model results accurate in 
reproducing fuel cell behaviour and validated for prediction, 
maintaining experimental uncertainty as prediction confidence 
interval. This implies that the phenomena neglected with the 
proposed assumptions have minor effects. Measures having 
lower uncertainty are needed for a further evaluation of model 
accuracy. Instead above 0.5 V the model results less accurate, 
implying that predictions in that range have an error higher than 
experimental uncertainty. 


4.3. Estimated quantities discussion 


The model estimations of methanol crossover fluxes are 
qualitatively and quantitatively similar to measures reported in 
literature, Fig. 13. The order of magnitude and the approxi- 
mately linear dependence on current density are consistent with 
experimental observations [7,21,22]. 


Current Density [A - cm~°] 


0 0.1 02 0.3 0.4 0.5 
Current Density [A - cm~7] 


Fig. 12. Residuals in function of current density. Experimental uncertainty in 
solid line. 
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Fig. 13. Methanol crossover in function of current density. Methanol concen- 
tration: 3.25 wt.%; varying anode flow rate and fuel cell temperature. 


In Fig. 14, examples of local current density in function of 
channel length are reported. At low voltage current density is 
homogeneous along the channel, instead at high voltage it is 
more intense at fuel inlet, near 0.5 A cm, and varies consid- 
erably along x. At 0.45 and 0.5 V, the transition from kinetics 
controlled region to transport controlled region is evident in the 
slope variation around, respectively, 500 and 150mm. In these 
cases, current production near channel outlet results very limited 
as expected. In Fig. 15, where local methanol concentrations are 
reported in logarithmic scale, the same transition is appreciable 
at a methanol concentration about 1 x 1076 mol cm. 

Void fraction increases rapidly along the channel; in the cases 
considered in Fig. 16 it reaches the value 0.5 at maximum within 
30 mm and values very close to 1 at channel outlet. 

In Figs. 17 and 18, total methanol flux through the diffusion 
layer, methanol flux in gas phase and the fraction associated to 
current production are presented. In the first example, methanol 
crossover, the difference between total flux and current flux, is 
very limited and present considerably in the first third of the 
channel. Different is the second case, where methanol crossover 
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Fig. 14. Local current density in function of channel length. Methanol con- 


centration: 3.25 wt.%; anode flow rate: 0.47 g min™!; fuel cell temperature: 
353K. 
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Fig. 15. Methanol concentration in the electrode in function of channel length. 
Methanol concentration: 3.25 wt.%; anode flow rate: 0.47 gmin™!; fuel cell 
temperature: 353 K. 
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Fig. 16. Void fraction in function of channel length. Methanol concentration: 
3.25 wt.%; anode flow rate: 0.47 g min~!; fuel cell temperature: 353 K. 
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Fig. 17. Methanol fluxes in diffusion layer in function of channel length. 
Potential difference: 0.5 V; methanol concentration: 3.25 wt.%; anode flow rate: 
0.47 g min™!; fuel cell temperature: 353 K. 
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Fig. 18. Methanol fluxes in diffusion layer in function of channel length. 
Potential difference: 0.5 V; methanol concentration: 6.5 wt.%; anode flow rate: 
1.08 g min~!; fuel cell temperature: 353 K. 


is always more than 30% of the total flow in the diffusion layer, 
reaching 50% at channel inlet. Regarding gas phase transport it 
is important to observe that it seems to produce the predominant 
contribution, excluding channel inlet. Moreover, its importance 
is evident also in the second case, despite the cell is fed with a 
very large methanol stoichiometry; in both cases more than 80% 
of the flux through the diffusion layer in the second half of the 
channel is due to gas phase transport. 

These observations confirm the important contribution that 
methanol transport in gas phase has in feeding the electrode and 
emphasise the possibility to optimise related aspects, in particu- 
lar diffusion layers, to reduce methanol crossover and to enhance 
fuel cell performance. 


5. Conclusions 


The analysis presented in this work produces the following 
conclusions: 


e A systematic experimental analysis of operating conditions 
influence on anode polarization is presented, characterising 
measurements in term of uncertainty and reproducibility. 

e Influence of operating condition history on performance is 
observed: a hysteresis effect on current measured values, 
decreasing and increasing voltage. Experimental results con- 
firm that it arises from methanol accumulation; consequently 
this phenomenon has to be considered in evaluating DMFC 
performances uncertainty and reproducibility of measure- 
ments. 

e Methanol electro-oxidation kinetics results weakly dependent 
on methanol concentration. 

e The developed model has been calibrated in an extensive oper- 
ating condition range, utilizing 48 current density measures, 
coming from eight polarization curves, and determining five 
fitting parameters. 

e The proposed model reproduces accurately DMFC anode 
polarization in the investigated range. The obtained val- 
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ues of fitting parameters and local quantities estimation are 
consistent with literature and the residuals between model 
estimations and experimental results are inferior to measure- 
ment uncertainty. 

e Both experimental results and model estimations confirm the 
important contribution of gas phase methanol transport in 
anode feeding. Moreover. they emphasise the possibility to 
optimise anode flow rate in order to improve DMFC perfor- 
mance and reduce methanol crossover. 
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